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Abstract 

In this paper we have considered a red-shift based linear correction in deriva- 
tive of action in the context of vacuum f(R) gravity. Here we have found 
out that the linear correction may describe the late time acceleration which 
is appeared by SNela with no need of dark energy Also we have calculated 
the asymptotic action for the desired correction. The value of all solutions 
may reduce to de' Sitter universe in the absence of correction term. 

Keywords: modified gravity, cosmology 
1. Introduction 

The recent data coming from the luminosity distance of SuperNovae la 
(SNela) Q, wide surveys on galaxies j^j and the anisotropy of cosmic mi- 
crowave background radiation |3[ suggest that the Universe is undergoing an 
accelerating expansion. Large Scale Structure formation Baryon Oscil- 
lations [Hf and Weak Lensing j^] also suggest such an accelerated expansion 
of the Universe. Actually, identifying the cause of this late time accelera- 
tion is one of the most challenging problems of modern cosmology. Several 
approaches being responsible for this expansion, have been proposed in the 
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literature. A positive cosmological constant can lead to accelerated expan- 
sion of the universe but it is plagued by the fine tuning problem Q. The 
cosmological constant may be interpreted either geometrically as modifying 
the left hand side of Einstein's equation or as a kinematic term on the right 
hand side with the equation of state parameter w — — 1. Another approach 
can further be generalized by considering a source term with an equation 
of state parameter w < —1/3. Such kinds of source terms have collectively 
come to be known as Dark Energy. Various scalar field models of dark energy 
have been considered in literature . All the dark energy based theories as- 
sume that the observed acceleration is the outcome of the action of a still 
unknown ingredient added to the cosmic pie. In terms of the Einstein equa- 
tions, G^ u = xT^iu, such models are simply modifying the right hand side 
including in the stress-energy tensor with something more than the usual 
matter and radiation components. 

On the other hand, as a radically different approach, one can also try to 
leave unchanged the source side, but rather than modifying the left hand side 
of Einstein field equations. In a sense, one is therefore interpreting cosmic 
acceleration as a first signal of the breakdown of the laws of physics as de- 
scribed by the standard General Relativity (GR). Extending GR, not simply 
given its positive results, opens the way to a large class of alternative theories 
of gravity ranging from extra dimensions [i| to non-minimally coupled scalar 



fields [lOj. In particular, we will be interested here in fourth order theories 



11] based on replacing the scalar curvature R in the Hilbert-Einstein action 



with a generic analytic function f(R) which should be reconstructed starting 
from data and physically motivated issues. Also referred to as f(R) gravity, 
these models have been shown to be able to fit both the cosmological data and 
Solar System constraints in several physically interesting cases [12]. These 
theories are also referred to as 'extended theories of gravity', since they nat- 
urally generalize General Relativity. It has been predicted that the universe 
might have been appeared from an inflationary phase in the past. It is also 
believed that the present universe is passing through a phase of the cosmic 
acceleration. Vacuum solutions of f(R) gravity theories are one of interesting 



subjects which are obtained for constant Ricci scalar [13j, [14], [15|, ll6[], while 



it is possible to derive non constant curvature scalar vacuum solutions. 

In this paper we would like to note that f(R) gravity theory is a powerful 
approach to describe dynamical behavior of the Universe via an unusual 
approach. Actually we consider vacuum solutions of f(R) gravity. But there 
is a difference with other vacuum solutions. This way we do not assume 
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constant scalar curvature to obtain vacuum solution. Vacuum solutions of 
modified f(R) gravity would like to explain the late time phase transition 
of cosmological parameters like deceleration parameter without the need for 
dark companion of the Universe, just by pure geometry. 

The pioneering works on reconstruction of modified action through inverse 



method are done by the authors of Ref. [13|, [14|, [15|, [16|j. They developed a 
general scheme for cosmological reconstruction of modified f(R) gravity in 
terms of e-folding (or red-shift) without using auxiliary scalar in intermedi- 
ate calculations. Using this method, it is possible to construct the specific 
modified gravity which contains any requested FRW cosmology. The number 
of f{R) gravity examples is used where the following background evolutions 
may be realized: AC DM epoch, deceleration with subsequent transition to 
effective phantom superacceleration leading to Big Rip singularity, deceler- 
ation with transition to transit phantom phase without future singularity, 
oscillating universe. It is important that all these cosmologies may be re- 
alized only by modified gravity without the use of any dark components. 
In this essay, we try to reconstruct an appropriate action for the modified 
gravity through the semi-inverse solution method. We do not assume any 
FRW cosmology to reconstruct its related f(R) action. Our starting point 
is some modification in deriving from a generic action which is depended on 
red- shift. 

In section II, we have a briefer review of modified field equations. In 
section III, we introduce our model and its results in field solutions. In 
section IV, we study the evolution of deceleration parameter under considered 
model and its related dark energy Equation of State (EoS). In section V, 
we calculate an approximated value for correction parameter which is in 
accordance with SNela data from observational constraints. In section VI, 
we try to reconstruct the original action which may produce our desired 
corrections and in section VII we examine the local tests for the obtained 
action. Section VIII is conclusion of this paper. 

2. Modified field equations 

The action of modified theory of gravity is given by 

S = J d i xV^g[^f(R) + L m ], (l) 

where L m is the matter action such as radiation, baryonic matter, dark mat- 
ter and so on which we do not consider them in field equation. In this essay, 
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we consider the flat Friedmann Robertson Walker, (FRW) background, so 
that the gravitational field equations for f(R) modified gravity are provided 
by the following form 

- 3 V + 3-Ar + £/ = 0, (2) 
a a 2 

[- + 2*]f - 2-Rf" - R 2 f - Rf" ~y = 0. (3) 

a a 2 a 2 

where the overdot denotes a derivative with respect to t, and the prime 
denotes a derivative with respect to R, a(t) is the scale factor and H = 
a(t)/a(t) is the Hubble parameter. Eliminating / between Eqs. (j2J) and ([31 
results: 

- 2[- - £) 2 }f + h -Rf" - Rf" - R 2 f" = 0. (4) 
a a a 

Which can be changed in the form of: 

2HF - HF + F = 0, (5) 

where F = df/dR. Eq. (JS]) is a second order differential equation of F with 
respect to time, in which both of F and H are undefined. The usual method 
to solve Eq. fl5]) is based on definition of f(R). 

Changing the variable of the above equation from t to a new variable 
like the number of e-folding, iV was done in 

HI 

Il5l . 1 1 61 ] . The variable N is 



related to the redshift, z by e~ N = 1+z. They solve the cosmological dynamic 
equation by definition of Hubble parameter as a function of N in a general 
form. Then they rewrite the equation by redefinition of the variable from 
the number of e-folding to the Ricci scalar and solve it with respect to the 
Ricci scalar. Thus they could demonstrate that modified f(R) gravity may 
describe the ACDM epoch without any need for introducing the effective 
cosmological constant, non-phantom matter and phantom matter. In this 
paper we would like to replace the variable of Eq. flSJ) by red-shift, z, directly 
to solve the dynamical equation for a specific red-shift depended action. 

Each redshift, z has an associated cosmic time t (the time when objects 
are observed with redshift z emitting light), so we can replace all the differ- 
entials with respect to t by z via: 

d dadz d 
dt dt da dz ' 

= _(1 + (6) 
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where we use 1 + z = a /a, and we consider a = 1, in the present time. 
Now, we can replace the variable of Eq. dSJ) from t to z by using Eq. 
and we obtain a first order differential equation for H 2 with respect to z as: 

±-H{zf = P(z)H{z)\ (7) 

where P(z) is a function with respect to z, which depends on the definition 
of F with respect to z as: 

_ 2(1 + z)(d 2 F/dz 2 )+A(dF/dz) 
[) 2F-(l + z){dF/dz) " [) 

Now we may solve Hubble parameter that depends on the definition of F(z). 
This may happen by analytical calculations or numerical approaches which 
is based on selection of F(z). Since there could be many choices to select the 
function of F(z), we have decided to add a linear correction term to General 
Relativistic limit of F(R) value, which is discussed in the following section. 

3. Model selection and its solution 

Since we do not know the determined function of F(z), would like to 
consider a linear correction of red-shift as 

F{z) = F (l + az), (9) 

in which for a = we have F(z) = Fq. In this case Eq. ((7j) has a solution 
such as H(z) 2 = Hq in which Hq is a constant. This case will reproduce de' 
Sitter solution of vacuum Universe that is expanded with a constant velocity. 
Since the value of H is independent of F we may fix it to Fq — 1 to recover 
GR or transfer its effect to gravitational constant by normalization. For 
a^Owe have our linear correction. In this case Eq. (J7]) has a solution such 
as 

H{zY = Hl{l-\a{l-z)]\ (10) 

where H 2 = H$/(l - a/2) 4 . It is clear that it may reduce to de' Sitter case 
in the absence of a. 
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4. Cosmic dynamics and EoS 



Deceleration parameter, q in cosmology is a dimensionless measure of the 
cosmic acceleration of the expansion of space in a FRW universe. It is defined 
by: 

da H . 

9 = -a>=-i-ip- <") 

here we change the variable of Eq. (lllj) from t to z then we will have evolution 
of deceleration parameter with respect to redshift as: 

1 + z dH(z) 2 , , 

= -l+ l -^P(z), (13) 

which only depends on P{z) for vacuum solution. Since we would like to 
study the first order correction of F(z) we put the modified model in Eq. 
ffTB"]) to obtain 

Q(Z) = "I + ^ -J-, ^T- 14 

2 — a{l — z) 

Transition point from deceleration to acceleration phase (q(z) = 0) obtains 
as 

z T = , (15) 

a 

which is a constraint for our correction parameter as0<a<2/3to have 
a transition point in recent positive red-shifts, but it should not be close to 
zero because this value will shoot the transition to high red-shifts. On the 
other hand as it is shown in Fig. (TjQ) for all the values of a < 2/3 universe 
is under acceleration. Also evolution of deceleration parameter from high 
red-shifts to the future is shown in Fig. (j2j). Another result of Eq. (I14p is 
for large red-shifts, q(z) — > 1 which is the value of deceleration parameter in 
radiation era. 

The equivalent form of dark energy Equation of State (EoS) which cor- 
responds to selected F(z) is obtained as 

1 + z dH(z) 2 . . 

= -1+ Ml + Z) (17) 
i + 3 [2_ a (i-^)]' 
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Figure 1: q(z = 0) with respect to a. 
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Figure 2: q(z = 0) with respect to a. 
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Figure 3: U)(z — 0) with respect to a. 



which is reduced to de' Sitter EoS that is equal to u — — 1 when a = 0. Fig. 
(j3J) shows the range of the predicted value of EoS with respect to acceptable 
values of a. Also evolution of EoS from high red-shifts to low red-shifts is 
shown in Fig. (j4]). On the other hand Eq. ( TT7T) shows that for red-shifts that 
are large enough, co(z) — > 1/3 which is EoS in radiation era, while we have 
solved modified Friedmann equations for vacuum. 



5. Supernovae constraint 

In the following we would like to constraint our model parameters in com- 
parison with observational data sets. Another main late time cosmological 
constraint which is tested for the above simple example is SNela distance 
module with respect to redshift constraint. In this comparison we used the 



Union2 data set [20| which provide 557 SNela specifications. According to 
the value of numerical results which was gotten in comparison with SNela 
data set from usual y 2 algorithm, we discover that for a = 0.418, h = 0.697 
where h is the reduced Hubble parameter and x 2 = 545.329, so the ratio of 
chi square error to the number of freedom is about 0.983. The result of this 
calculation is shown in Fig. (JSJ). 



6. Reconstruction of f(R) action 

Specification of modified action is one of interesting subjects in the story 
of f(R) gravity. Here we may reconstruct the action via a semi inverse 
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Figure 4: Evolution of u(z) for KCDM (solid line) with il m = 0.3 and Oa = 0.7, and for 
our model with correction parameters a = 0.5 (dotted line), a = 0.4 (dashed line) and 
a = 0.3 (dashed-dotted line). 
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Figure 6: Evolution of R/Rq with respect to z for three values of correction parameter. 
The crossing points are not unique. 



solution approach. The Ricci scalar is defined as 

R = Q(H + 2H 2 ), 



(18) 



and if we change the variable of this equation from time to red-shift we have 

dH(z)\ 



R(z) = 3[4H(z) 2 - (1 + z) 

The Ricci scalar for the Hubble parameter of Eq. 
(DEI) as 



dz 



(19) 



R{z) = R a 



l-l a (l-z) 



TUT) is obtained from Eq. 



(20) 



where 



Rn 



12(1 



a)Hl 



(l-a/2) 4 ' (21) 

which will be reduced to another constant in a = and we may name it 
the present value of Ricci scalar Ro = 12Hq. Changes of modified Ricci 
scalar is shown in Fig. ([6]). If we suppose the value of cosmological constant, 
A = 3Hq then we have a simplified relation between the present value of 
curvature scalar and cosmological constant as Rq = 4A. We may use these 
equalities to find behavior of the modified solutions in the constraint regions. 



10 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 

R 

Figure 7: Change of f{R) with respect to R for three values of correction parameter, a = 
in dashed line, a — 0.3 in dot-dashed line and a — 0.6 in dotted line. Solid line shows 
variation of GR action with cosmological constant which is de' Sitter action. 

One of the simple solutions to achieve the desired action is eliminating of 
z between F(z) and R(z) which leads us to reach the equation 

F( J R) = 2(A)i/3 + a _ 1 . (22 ) 

It is clear that, the derivative of action is different from GR plus cosmological 
constant or de' Sitter model. But in the case of a = or R a = R , we would 
like to recover de' Sitter solution. Then by using Eq. ( 120]) we obtain R = Rq 
and then F(R) = 1 which is equal to the GR or de'Sitter model. Because 
GR action differs from de'Sitter action only in a constant term. 

Now we can calculate f(R) by integrating of Eq. ( J22l with respect to R 

as 

/(^) = ^(^) 1/3 + «-l] + C, (23) 

where C is constant of integration. Since we would like to recover de' Sitter 
model in the absence of correction term, a we must find the true constraint on 
the constant of integration, we can determine the constraint by comparing 
the value of de' Sitter action with the action of Eq. ( |23l) when a = or 
R = Ro. The standard form of de' Sitter action is fd(R) = R — 2A. Since the 
scalar curvature of de' Sitter universe is a constant, then we can determine 
the value of fd{R) in the presence of constant curvature. If we solve the 
standard algebraic equation of f(R) gravity as 2f d (R) — RF d (R) = we can 
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find Ro = 4A. Then the value of de' Sitter action in the case of R = Rq is 
f d (R = 4A) = 2A which is equal to R /2. 

For the case of our obtained action from Eq. f )23|) . we have f(R = Ro) = 
Ro/2 + C when a = 0. Since we would like to have the same solution for de' 
Sitter model and our obtained action in present of R = Rq we force constant 
of integration, C to be equal to zero. 

Therefore the obtained action in its constrained form would be 

/(^) = ^(^) 1/3 + «-l], (24) 

which is different from de' Sitter action even in the limit of a = 0, but it 
behaves like fd{R) = R — 2A in the region of R ~ Ro as shown in Fig. (JZJ). 
In the following of this section we discuss one special case of parameter a to 
see how the reconstructed action is related to the de' Sitter action. 

6.1. a < 1 

In this case the constant R a is reduced to R a ~ Ro(l + a). There are 
two subcases in this limit. The first one is low red-shift approximation. In 
this case Ricci scalar in Eq. ( 120]) is reduced to R(z) ~ Ro(l — at/2) = const., 
which is as same as the case of constant expansion rate or de' Sitter universe, 
but reconstruction of F(R) is impossible via our approach. The second one is 
high red-shift approximation. In this case Ricci scalar is reduced to R(z) ~ 
Ro(l + \ctz). Then the derivative of action is obtained in the form of 

F(J2) = ±(l+2-£). (25) 

Integration of the above equation with respect to R obtains the desired action 

as 

W) = kn + $-) + c, (26) 

where C is the constant of integration which plays the role of effective cosmo- 
logical constant. The factor of 1/3 may transform the gravitational constant 
and translate Gn to G e ff. In this case the form of the above action is reduced 
to de' Sitter action when R — > and effective cosmological constant defines 
as A e ff = 3C. 
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7. Local tests 



Following 17| here we introduce the auxiliary field A to rewrite geometric 



part of the action ([I]) in the following form: 

s = yJ d * x ^~ 9 {f{A) {R ~ A) + f{A)} ■ (27) 

By the variation over A, one obtains A = R. Substituting A = R into the 
action (|27|) . one can reproduce the action in (pQ). Furthermore, we rescale the 
metric in the following way (conformal transformation): 

g^^e'g^, a = -\nf(A). (28) 

Hence, the Einstein frame action is obtained: 

S E = ±-J d 4 x^g (r - ^g p °d p ad a a - V{a)\ , (29) 

A f{A) 



V{a) = e ^(e- CT )-e 2< 7(<?(e- CT )) 



f'(A) f'(Af 



Here g (e _CT ) is given by solving the equation a = In f'(A) as A = g(e~ a ). 
Due to the scale transformation (|28|) . a coupling of the scalar field o with 
usual matter appears. The mass of a is given by 



2 _ ld 2 V(a) 1 f A 4/(A) 



m ^-2 ^ 2 \f (A) (f(A)) 2 + f"(A) ) ■ (30) 

Unless m a is very large, there appears large correction to the Newton law. 
Naively, one expects the order of the mass m a to be that of the Hubble rate, 
that is, vn a ~ H ~ 10~ 33 eV, which is very light and the correction could be 



very large, which is claimed in [17 



We should note, however, that the mass m a depends on the detailed form 
of F(R) in general }18j . Moreover, the mass m a depends on the curvature. 
The curvature on the earth -R ea rth is much larger than the average curvature 
-R so i ar in the solar system and -R so iar is also much larger than the average 
curvature in the unverse, whose order is given by the square of the Hubble 
rate H 2 , that is, R e£LT th -Rsoiar 3> H 2 . Then if the mass becomes large when 
the curvature is large, the correction to the Newton law could be small. Such 
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mechanism is called the Chameleon mechanism and proposed for the scalar- 
tensor theory in 19|. In the case of action §M§, the mass m a is given by 



ml ~ \r 1 1*A*I\ (31) 

where R a is defined in Eq. f l2Tj) . The order of R a is about R a ~ 12Hq ~ 
10 _65 eV 2 for the value of < a < 2/3. Then in solar system, where R ~ 
10 _61 eV 2 , the mass is given by m 2 ~ 10~ 31 eV and in the air on the earth, 
where R ~ 10~ 50 eV 2 , m a ~ 10 _28 eV. The order of the radius of the earth 
is 10 7 m ~ (10~ 14 eV) , m a ~ 10~ 13 eV. Therefore the scalar field a is very 
light and the correction to the Newton law is observable. Here we should 
note that the obtained action ( 124"|) strongly depends on the order of red-shift 
based correction terms, strongly. Then one can consider local tests under 
influence of higher order corrections in the derivative of action. 

8. Conclusions 

In this paper we have examined a first order red-shift based correction 
on derivative of general action with respect to Ricci scalar as a starting 
point of modified gravity in f(R) theory. Also we have shown that this 
correction may operate as an alternative for dark energy via its ability in 
according with distance modulus of Union2 data set of SNela. There is an 
interesting behavior in the limit of low curvature regions for the obtained 
modified action. As it is clear in the region of R — > of Fig. (J7|), the 
modified action tends to the value of scalar curvature, f(R) — > R for all 
of values of correction term, which is pure General Relativity action and 
in the limit of R — >■ Rq it will be close to f(R) — > R — 2 A in the case of 
a = 0. However the mass of equivalent scalar field of the obtained action 
is not heavy enough to evade local tests of the theory, the behavior of other 
spacetime solutions such as spherically symmetric space is considerable. As 
it is shown in this paper, solution of modified equation in the absence of 
baryonic matter and cosmic cold dark matter could get better results in 
comparison with AC DM in the case of expansion rate of universe. In this 
manner the f(R) may play as an alternative for cosmic dark matter. But 
f(R) theories can also play a major role at astrophysical scales. In fact, 
modifying the gravitational Lagrangian affects the gravitational potential 
in the low energy limit. A corrected gravitational potential could offer the 
possibility to fit galaxy rotation curves without any need of huge amounts 
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of dark matter, which is considerable for the obtained action (I24p . On the 
other hand behavior of high red-shifts or early time cosmology of the action 
is considerable, because the power law terms may be considered to explain 
inflationary behavior of early universe. Of course the main correction of 
this approach comes from selection of F(z) = 1 + az which is a toy model 
that could to explain late time inflationary behavior of the universe without 
dark energy. This approach is based on general behavior scalar curvature, 
not as a constant. Also it will explain radiation dominated era values of 
parameters such as deceleration parameter and EoS which is in accordance 
with a universe including radiation density, while the modified field equations 
have been solved for empty universe. Our main scope from this paper is 
introducing an approach that may proposed as an alternative for dark energy. 
So we can demonstrate if the obtained action is viable for other scales of 
universe or not? Finally we propose higher order red-shift based corrections 
to have more accurate behavior of modified solutions which is in progress. 
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